Electron-Electron Interactions in Graphene 
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We discuss the validity (or not) of the ring-diagram approximation (i.e. RPA) in the calculation 
of graphene self-energy in the weak-coupling (r a <C 1) limit, showing that RPA is a controlled and 
valid approximation for extrinsic graphene where the Fermi level is away from the Dirac point. 
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We point out in this Comment that a recent Letter 
[l[ by Mishchenko is both incorrect and misleading in 
asserting that the RPA self-energy is not a controlled ap- 
proximation for graphene in spite of graphene being a 
weak-coupling two-dimensional (2D) system in the sense 
that the interaction coupling constant r s — e 2 /(Hkv), 
where k, v are respectively the background lattice di- 
electric constant (k ~ 2.5 for graphene on SiO^) and 
the graphene band velocity (v ~ 10 6 m/s), is small (i.e. 
r s < 1) in graphene. We emphasize here that for extrinsic 
graphene (i.e. gated or doped graphene with a free car- 
rier density-induced chemical potential or Fermi level Ef 
in the conduction or valence band away from the Dirac 
point) the RPA is an excellent and controlled approxi- 
mation which is asymptotically exact in the r s — > limit 
precisely as it is in a 2D or 3D electron gas in a jellium 
background @, I n intrinsic graphene, however, where 
the system is a zero-gap intrinsic semiconductor with the 
chemical potential precisely at the Dirac point the 
RPA self-energy calculation, for any arbitrary coupling- 
constant value (including r s -C 1), is known [H, |(| to 
lead to the failure of the Fermi liquid theory with the in- 
teracting system becoming a marginal Fermi liquid with 
a logarithmically divergent quasiparticle rcnormalization 
factor, and as such all perturbative many-body approxi- 
mations become suspect, including the RPA. The sweep- 
ing statement made in Ref. dj about the lack of validity 
of RPA in calculating the graphene self-energy, even in 
the weak-coupling r s <C 1 regime, is thus incorrect for ex- 
trinsic graphene and trivial for intrinsic graphene. The 
conceptual confusion in Ref. [1] arises from an inabil- 
ity to distinguish between the large-momenta ultraviolet 
divergence and the low-momenta infrared divergence as 
described below. 

In Fig. 1 we show the leading-order many-body per- 
turbative Feynman diagrams to the graphene self-energy 
with the bare electron-electron interaction being the 
long-range Coulomb interaction V c (q) — 2ire 2 / (nq). The 
first-order (in V c ) diagrams, given by Figs. 1(a) and 1(b), 
are simple: The Hartree tadpole diagram (Fig. 1(a)) van- 
ishes by virtue of the charge neutrality due to the lattice 
background and the Fock exchange diagram, Fig. 1(b), 
has been calculated and discussed in details in the lit- 
erature p|. Of the three diagrams in the second order, 



FIG. 1: First-order [l(a)-l(b)] and second-order [l(c)-l(e)] 
self-energy diagrams. Fig. 1(f) shows the self-energy in the 
infinite ring-diagram approximation (i.e. RPA). The solid 
line denotes the electron Green function and the wavy line 
the bare Coulomb interaction. 



only Fig. 1(d) and (e) are non-zero since the contribution 
of Fig. 1(c) can be trivially subsumed into the exchange 
self-energy diagram of Fig. 1 (b) . For the usual parabolic 
dispersion jellium background electrongas problem dis- 
cussed extensively in the literature @, [J, the long wave- 
length infrared (i.e. q — ► 0) singularity of the Coulomb 
interaction dominates the energetics in the r s <C 1 limit, 
and it is well-known that the ring diagram of Fig. 1(e) 
has a logarithmic divergence which is fixed by summing 
the subset of higher-order diagrams shown in Fig. 1(f) 
which is the infinite series of ring diagrams. This is 
the RPA technique which is equivalent to keeping the 
leading-order term in the renormalized expansion using 
the effective dynamically screened Coulomb interaction 
u (l) = Vc{q)/e(q,u) where e(q,u>) = 1 - V c (q)U(q, u), 
with II being the bare electron-hole ring or bubble dia- 
gram or equivalently the non-interacting irreducible po- 
larizability function. Note that e arises from the infi- 
nite geometric series of keeping all the bubble diagrams 
[Fig. 1(f)]. The RPA re-summation technique completely 
regularizes the Coulomb infrared divergence problem and 
is therefore asymptotically exact in the r s — > limit. 

Graphene, in addition to having the Coulomb infrared 
divergence, also has an ultraviolet divergence (i.e. q — > 
od) arising from its peculiar chiral linear band dispersion. 
This ultraviolet divergence is usually regularized by in- 
troducing a physical momentum cut-off q c ~ a -1 where 



2 



a is the graphene lattice constant. This typically intro- 
duces logarithmic terms in the cut-off momentum. This 
ultraviolet divergence exists for all the self-energy dia- 
grams for graphene, and there is nothing special about 
the ring diagrams as far as the ultraviolet divergence is 
concerned in contrast to the infrared divergence where 
the ring diagrams dominate in the r s <C 1 limit. For the 
ultraviolet divergence itself, the ring diagram [Fig. 1(e)] 
is not in any sense more important (or, parametrically 
larger) than the vertex diagram of Fig. 1(d) for any value 
of r s . 

The mathematical point for graphene, erroneously ig- 
nored in Ref . [H , is that while intrinsic graphene has only 
the ultraviolet divergence (where the ring diagrams and 
RPA have no special significance) , extrinsic graphene has 
both the infrared and the ultraviolet divergence, mak- 
ing RPA a well-controlled approximation in the r s <C 1 
regime. Thus, RPA is a perfectly meaningful approxima- 
tion for extrinsic (i.e. doped or gated) graphene where 
much of the experimental work is being done. 

To demonstrate the subtle mathematical point, we 
consider the potentially singular contribution E(fc,w) 
to the graphene self-energy coming from the ring di- 
agram of Fig. 1(e), finding it to go as T,(k,oj) = 
r 2 J °° dql(q,k,uj)/q where, for intrinsic graphene, 
I(q, k, lj) ~ q as q — > and ~ constant as q — * oc. Thus 
intrinsic graphene has only the ultraviolet, but not the 
infrared, divergence, and E(fc) ~ hi(q c /k), implying that 
ring diagrams or RPA have no special significance for in- 
trinsic graphene. For extrinsic graphene, however, I(q) ~ 
constant for both q — > and q — > oo, and therefore, the 
infrared (as well as the ultraviolet) divergence is present, 
making RPA a meaningful well-controlled approximation 
0,1, as in ordinary electron gas systems in the weak- 
coupling regime. In extrinsic graphene, which is the ex- 
perimentally relevant graphene system with free carriers, 
RPA regularizes the usual infrared divergence associated 
with the long-range Coulomb interaction whereas RPA 
is not meaningful in intrinsic graphene which does not 
manifest any infrared Coulomb divergence. More details 
are given in the appendix. 

We point out that in graphene, even extrinsic graphene 
where the carrier density can be varied, the coupling 
parameter r s — e 2 /(hnv) does not depend on the den- 
sity, in sharp contrast to the usual 2D (r s ~ n -1 / 2 ) 
and 3D (r s ~ n -1 / 3 ) parabolic dispersion electron sys- 
tems. Thus, in graphene r s ranges between 2.4 (for free 
standing graphene with k — 1) and — e.g. r s ~ 0.9 
for graphene on Si02, r s ~ 0.15 for graphene on Hf02- 
Therefore, in some sense RPA is more accurate in ex- 
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trinsic graphene than in 3D metals (r s w 3 — 6) since 
extrinsic graphene (with r s w 0.15 — 2.4) better satisfies 
the weak-coupling (r s < 1) condition necessary for the 
quantitative validity of RPA. 

Finally, we also mention that the main result reported 
in Ref. [l|, namely that the intrinsic conductivity in 
graphene is renormalized by electron-electron interac- 
tions, is manifestly incorrect as it violates the well-known 
Kohn theorem. The long wavelength conductivity cannot 
be modified by electron-electron interaction in a trans- 
lationally invariant system by virtue of the separation 
of center of mass and relative coordinates, and the ver- 
tex correction and the self-energy diagrams must cancel 
each other in each order, as shown recently explicitly for 
graphene in Ref. [§]. 

APPENDIX 

The kinetic energy of graphene for 2D wave vector k is 
given by e s k = su|k|, (we use H = 1 throughout this pa- 
per) where s = ±1 indicate the conduction (+1) and va- 
lence (—1) bands, respectively, and v is the Fermi veloc- 
ity. The corresponding density of states (DOS) is given 
by D(e) = g\s\/{2Trv 2 ), where g = g s g v with g s = 2, 
g v = 2 the spin and valley degeneracies, respectively. 
The Fermi momentum kp and the Fermi energy ep of 2D 
graphene are given by kp = (4nn/ g s g v Y^ 2 and £p = vkp 
where n is the 2D carrier (electron or hole) density. For 
intrinsic (extrinsic) graphene, n, ep = (t^)O with the 
Dirac point taken to be the energy zero. 

The exchange energy graph is the first-order diagram 
of interaction shown as in Fig. 1(b). It is the only contri- 
bution with one Coulomb line. The exchange self-energy 
is given by 

Ef (k) = n F (& +qs ') VM) F$ (k, k + q), (1) 

where £ ks = e s k-£F, n F (£ ks ) = 6(kp-sek) is the Fermi 
function at T = 0, F^) (k, k + q) is the band overlap 
matrix element given by 

F s ( s V(k,k') = i(l + SS 'cos0 kk O, (2) 

where k k' is the angle between k and k'. 

There are three diagrams in the second-order self- 
energy diagrams with two Coulomb lines. Fig. 1(d) can 
be expressed as 



EW(k,tA n ) = -l E E E K(q)K(q0^ 2 L S 3( k > k i> k 2,k 3 )G^^ (3) 

si ,S2 1S3 q-q' iq n ,iq' n 
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where (3 = l//sgT, s, Sj = ±1 denote the band indices, G°(k, ik n ) = l/(ik n — £ks) is the unperturbed Green function, 
V c (q) = 2ne 2 / nq is the bare Coulomb potential (with background dielectric constant n). And the function is the 
band overlap matrix element, given by 



Fi i 2 ) 2S3S4 (k 1 ,k 2 ,k 3 ,k 4 ) = i 



I 4 

1 + S1S2S3S4 + -Z ^2 SiS i C0S ^i.kj 



After summing over frequencies and a standard procedure of analytical continuation, we have 

Vc(q)W) 



SW(k,c) 



E E 



-Fi a 2 1 ) a2a3 (k,k 1 ,k 2 ,k3) 



x {nF(^k 3S3 ) [«F(Ck lS i) - n F (& 2S2 )} + n F (& 2S2 ) [1 - n F (£k lSl )]} 

= E E fi ™( k >^q><i')- 



5 1 ,s 2 ,S3 q.q' 



(4) 



(5) 



For intrinsic graphene nF(£k+) 

Thus, only the following terms in are nonzero for 
intrinsic graphene: i_ and R 



and n F (£k-) = 1. Eq. (5): R { +1 ++ , R { f ++ 
Here R ( f +++ , <1+ 
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do not show any sin- 
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As 



00, 



these functions behave as \/q. Therefore the self-energy 
contributions by these terms are logarithmically diver- 
gent. After introducing a physical cutoff (q c ) we have 
Eg^ oc Cdhi(q c /k), where d is a regular function of 

k,oj. As q — > the integrands R + ^ |_ and R^ | are 

finite and introduce no singular behaviors in E s d \ 

For extrinsic graphene nF(£k+) = 9(kp — |k|) and 
^F(^k-) = 1- Thus, we have five non-zero terms in 



and Rf 



and q — ► 00 limits. Again 



gular behavior in both q 

__i show singular behavior, ocl/q, only 

Thus we can control this singularity by intro- 



++- 
as q — > 00. 

ducing a physical cutoff. In summary of E s d \ we have 
an ultraviolet divergence in the integrand both for intrin- 
sic and extrinsic graphene. However, we can remove this 
divergence by introducing a cutoff q c <~ 1/a. 

Fig. 1(e) can be written as 



" si q jg n 

where n(q,ig n ) is the polarizability (bare bubble diagram) of graphene, given by 

n (q , ig„) = fl £E M^-M^) ^ (k , k + q) . (7) 

si,«2 k 4 <M ' Sksi ?k+qs 2 

We can express Eq. (6) as E^ (k, w) = E s lne (k, u) + EP olc (k, lu) after a standard procedure of analytical continuation: 
ES-OC*) = -EE / ^[^(q)] 2 e)(k,k + q) . , (8) 

SP° lo (k, u>) = £ ^ - W') - 0(-&+q.')] [K(q)] 2 F«(k, k + q)H(q, £ k+qSl - w). (9) 

Siq 
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Here EP olc does not show any singular behavior and is S s lne as 
a well-defined function for all k and w. We can rewrite 



£f e (k,u,) = -r s 2 - jf -fl(q,k,u). (10) 
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For intrinsic graphene we have I(q) ~ q as q — > 0, and 
I(<7) ~ constant as q — > oo. Thus, we have only the ultra- 
violet divergence in the integrand, but no infrared diver- 
gence. Again, si 6 - 1 oc ln(g c /fc). For extrinsic graphene we 
have a different behavior as q — > 0, i.e. 7(g) ~ constant 
as q^O. Therefore, we have both infrared and ultravi- 
olet divergences in the integrand. This singular behavior 
arises from the 1/q nature of the Coulomb interaction. 
We cannot control this divergence by introducing a cut- 
off, and the q — ► divergence must be regularized by the 
usual infinite ring-diagram approximation as in the usual 
RPA [i.e. Fig. 1(f)]. 

In conclusion, we have the ultraviolet divergence con- 
trolled by a cutoff for all diagrams but only Fig. 1(e) 
(RPA-type diagram) shows the infrared divergence for 
extrinsic graphene, which has to be controlled by the 
RPA re-summation given in Fig. 1(f). 



[1] E.G. Mishchenko, Phys. Rev. Lett. 98, 216801 (2007). 

[2] A.L. Fetter and J.D. Walecka, Quantum Theory of Many- 
Body Systems (McGraw-Hill, New York, 1971). 

[3] G.D. Mahan, Many-Particle Physics (Plenum Press, New 
York, 2000). 

[4] E.H. Hwang, Ben Yu-Kuang Hu, and S. Das Sarma, 

|arXiv:cond-mat/0703499V l. 
[5] S. Das Sarma, E.H. Hwang, and Wang-Kong Tse, Phys. 

Rev. B 75, 121406(R) (2007). 
[6] J. Gonzlez, F. Guinea, and M.A.H. Vozmediano, Phys. 

Rev. Lett. 77, 3589 (1996). 

[7] E.H. Hwang et al, |cond-mat/0 612345 Y. Barlas et al, 

Phys. Rev. Lett 98, 236601 (200 7). ~ 
[8] M. Poli ni et ffl/.,larXiv:0707.4230l E.H. Hwang and S. Das 

Sarma, larXiv:0708.1~133l 

[9] D.E. Sheeny and J. Schmalian, rarXiv:0707.2945V 2. 



